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Abstract 

By using the inverse spectral transform, the SRS equations are solved 
and the explicit output data is given for arbitrary laser pump and Stokes 
seed profiles injected on a vacuum of optical phonons. For long dura- 
tion laser pulses, this solution is modified such as to take into account the 
damping rate of the optical phonon wave. This model is used to interprete 
the experiments of Driihl, Wenzel and Carlsten (Phys. Rev. Lett., (1983) 
51, 1171), in particular the creation of a spike of (anomalous) pump radi- 
ation. The related nonlinear Fourier spectrum does not contain discrete 
eigenvalue, hence this Raman spike is not a soliton. 
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Stimulated Raman scattering (SRS) is one of the most studied three-wave 
interaction processes in nonlinear optics not only because it retains all the in- 
gredients of any other stimulated process, but also because it has revealed many 
striking, and sometimes unexplained phenomena. The theory of SRS has been 
developped on a semi-classical basis, for instance by Shen and Bloembergen jjj, 
Wang H or Carman et al Q| by assuming a permanent pump intensity pro- 
file. But when "depletion of the laser power implies that the laser field may not 
be treated as a fixed constant parameter" B, Stokes generation and amplifica- 
tion induces pump depletion and it is a serious obstacle to the propagation of 
high intensity laser pulses in a Raman-active medium. However it can also be 
thought of as a means to study (experimentaly and theoreticaly) the fundamen- 
tal properties of matter and radiation, and indeed Raymer and Mostowski |^| 
predict large (macroscopic) fluctuations of the Stokes pulse energy that are rem- 
iniscent of the small (quantum) fluctuations of the material dynamical variable 
(the polarization density state variable). These predictions were then cheked 
on measurements of the statistical distribution of Stokes pulses in pressured Hi 
gas by Wamsley and Raymer || and also by Fabricius, Natterman and von der 
Linde §. 

At the same time, experiments of Driihl, Wenzel and Carlsten [Q for long 
duration pump pulses (order of 100ns) revealed that "ocasionally the pump 
depletion is anomalously reversed for a short time interval generating a spike 
of pump radiation". They call it a soliton refering to the soliton-solution of 
the undamped SRS equations given by Chu and Scott |8| . The lesson of these 
results is the spectacular fit of the experimental data with numerical simulations 
of the SRS equations, indicating that the model is quite adequate. Then the 
fundamental result is the discovery that the spike in the pump-depletion occurs 
when a 7r-phase shift is introduced in the Stokes seed. 

In the absence of phase flip, the spike formation is still obseved (on real 
experiments, not on simulations) but with much lower probability. This led 
Englund and Bowden J^] to propose that the spike is the macroscopic man- 
isfestation of quantum phase fluctuations of the Stokes wave. The subsequent 
experiments of MacPherson, Swanson and Carlsten revealed that the anomalous 
pump radiation spike (what they call a Raman soliton) occurs in about 10.1% 
of the shots j|o|. Englund and Bowden |ll| developed a complete theoretical 
basis of such Raman spikes generated by quantum fluctuations of the initial 
Stokes vacuum, and they obtained a reasonable qualitative agreement with the 
previous experiments (they found a spontaneously generated spike in 13.6% of 
the shots). 

The observed spikes of pump radiation acquire then importance also for 
fundamental studies and they are always refered to as being solitons |l2| . But 
it was already remarked on the first series of numerical simulations that the 
spike narrowing (as propagation distance is increased) indicates that they are 
"not solitary waves in the strict sense" M. In particular the non-zero velocity of 
the Raman soliton implies, as shown by Menyuk, that "any solitonlike structures 
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are subluminous and will ultimately disappear at the back end of the pulse" fll3[ . 

Consequently the problem of the theoretical interpretation of the experi- 
ments of Driihl, Wenzel and Carlsten |7| is still open and we solve it here in terms 
of the inverse spectral transform theory (1ST) extended to arbitrary boundary 
values [jl5). We obatain the explicit global solution (output laser pulse) which 
maps perfectly the numerical simulation of SRS (as performed by MacPherson, 
Carlsten and Druhl 0). 

Such an explicit analytic formula for the output (eqs. (0) and (|^) below) is 
important for physics in many aspects: 

1 - It allows to understand why the Stokes phase flip and the finiteness of the 
dephasing time (long pulses) are both essential to the spike formation, and to 
discover the precise nonlinear mechanism generating the Raman spike. 

2 - It provides a powerfull tool to analyse the experiments. Indeed, having the 
digitalized input, our formula readily gives the predicted output and a compar- 
ison with the data of JjJ will be published later with other details [H . 

3 - It unveils the nature of the Raman spike as its nonlinear Fourier spectrum 
does not contain isolated points (bound states) but consists only in the con- 
timuum (radiation). Apart form the fact that the spike is not a soliton, the 
important consequence of this is that it survives long propagation distances, 
which is important for applications. 

The last point to mention is the fact that our solution, although being 
obatincd from the infinite line case, is quite close to the finite line case solution 
(obtained through numerical simulations). This is easily be seen by comparing 
for instance our figure 2 to the figure 2 of |l7|. The reason for this is that the 
input data are of finite duration and the interaction is local (for any fixed time) . 
We show in that indeed the generated material excitation is localized in a 
small region (of length comparable with the duration of the amplified Stokes 
pulse) . 

The model of SRS can be taken for instance from and reads, if we neglect 
the ground state depletion (i?3 = — N): 



The spatial variable ( lies in [0, L] , where L is the total beam path in the Raman 
cell, and the retarded time r = t — £/c is positive. Al and As are the slow 
envelopes of the laser pump (frequency u>l) and of the Stokes emission (fre- 
quency to s) which stimulate the material dynamical variable R (optical phonon, 
frequency up — lul — u>s)- K^s is the complex coupling constant, T<i the mean 
collisional dephasing time and N is a scaled density (= pAL with p the density 
of Raman active molecules and A the effective crossectional area of the pump 
beam). The initial-boundary value problem associated to the system (Q) is the 
following: 



d c A L = K LS RA S , d c A s = -K* LS R*A L 
8 T R + RL/(cT 2 ) = -K LS NA L A* S . 



(1) 



i?(C,0) = 0, A l (0,t)=A lo (t), A s (0,t) = A so (t), 



(2) 
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where, to reproduce the experiments, Alo is a gaussian, Aso is a fraction of 
Alo with possibly a change of sign (phase flip) somewhere. The problem is to 
determine the output quantities Al(L,t) and As{L,t). 

Our model results from (|l|) by considering first an infinite line (L — > oo) and 
taking into account the mismatch wave number of value 2k = kp+ks — k^, which 
results for instance from the Doppler effect due to molecular thermal motions. 
The resulting system is ([l]) but with Ag replaced with As exp[— 2ikQ. To take 
into account the contributions of all values of k, we introduce the distribution 
g(k) (centered in k — 0) of the relative coupling intensities. The resulting model 
equations are 

d x ai = qa 2 , d x a 2 — 2ika 2 = —qa\ 

dtq + iq = j dkg(k)aia 2 - (3) 

Here and in the following, an integral with no specified boundaries stands for 
(— oo,+oo). Hcrcabove we have made the following change of variables and 
scalings: 

x = -C, t = T, q = -K LS R, j = L/(dT 2 ), (4) 

(fc = 0): ai =A L /Ao, a 2 = A s e 2lkx /A , (5) 

where Aq = Mclx{Al(Q, t)}. The distribution g(k) is actually related to the 
inhomogeneous broadening and can be normalized to the coupling constant by 
setting J dkg(k) = K1 s N\Aq\ 2 . The related initial-boundary value problem 
corresponding to (||) reads here (note the sign in (0)) 

q(x,0)=0, ai(fc,+oo,f) = h{k,t), a 2 (k, +oo, t) = I 2 (k, t)e 2lkx (6) 

where for k = 0, h(t) = A lq (t)/A and I 2 (t) = A S o(t)/A . The problem to be 
solved is now to compute the output data a\[k, —oo, t) and a 2 [k, — oo, t). 

Note that, although t represents the retarded time, the initial value problem 
is physically meaningful because, the medium being initially in the ground state, 
we have set q(x, 0) = 0. Another important remark is that the model (||) maps 
onto (jl]) in the limit when g(k) becomes the Dirac distribution 5(k). However, 
on the infinite line, this is a singular limit ||l5f (in short it is not compatible 
with q(x,t) — > as x — * ±oo), and hence g(k) can be as sharp as we want but 
never a true delta function (which is physically quite reasonable). 

It has been shown in |lq| that the above initial-boundary value problem is 
solvable for 7 = and we shall use here directly these results. The solution is 
given by 

ai(fc,-oo,i) = h/P + I 2 a/0, (7) 

a 2 {k, -00, t)e- 2lkx =I 2 /p-ha/0. (8) 

The coefficients a(k,t) and (3(k,t) (the spectral data) can be computed ex- 
plicitely from and read for 7 = 

a(k, t) = -irg(k)e^ k ' t} f dt' I 1 (k,t')I 2 (k,t / )e-^ t "> (9) 
Jq 
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<t>(k,t) = / d£ 



lng(k)U(k,t')-~l^g(\)U(\,t') 



(10) 



U(k,t) = \h(k,t)\ 2 -\I 2 (k,t)\ 2 (11) 

p(k, t) = VTTRW e ie[k ' t] (12) 
m t) = --L I Y^k log(1 + |a(A ' <)|2) ' (13) 

where the slashed integral denotes the Cauchy principal value. Analogous for- 
mulae can be found in M ] but in the context of resonant interaction of light 
with a two-level medium (and application to superfluorescence) for which the 
boundary value problem notably differs. The above result gives the exact solu- 
tion to SRS for short pump pulses (for which 7 ~ 0). We report the discussion 
of this case to forthcoming paper and consider now the case of long pump pulses 
for which the pump depletion can be anomalously reversed . 

Our main argument is that both pump depletion and spike formation are de- 
scribed by the above solution of the boundary-value problem (^) for the system 
(||). Indeed, considering (Q) and (|l^), we remark that 

a — > 00 => |/3| — » 00 => |ai(-oo)| — » |7 2 |, (14) 

I /3 1 > 1 => |ai(-oo)| -> (15) 

Then pump depletion will occur in the time region where a(k, t) is large (the 
pump output will be of the order of the Stokes input I2), and pump radiation 
will occur in the time region where a(k, t) is close to zero (the pump output will 
be of the order of the pump input 

In the case of long duration pump pulses, the effect of the dephasing time 
is included in our model by assuming instead of the evolution (||) the following 
one 

a{k, t) = - vg (k)e^ k ^ / dt'I 1 (k,t')I 2 (k,t')e-' l > {k < t ' ) e' y( - t '- t l (16) 



(i 



The added exponential factor hereabove is justified by considering the linear 
limit in the spectral transform context. To that end it is convenient to write 
down the solution of (||) given by the spectral transform method [H : 

at =hh+he 2lkx h, (17) 

a 2 = -hh+he 2lkx h, (18) 

q = -^Jdkae~ 2lkx f 1 , (19) 

where f\ and f 2 are the solution of 

/ 1 (fc,x,i) = l + J- f J* a{\,t)e 2 ^f 2 {\x,t), 
Ziir J a + iv — k 
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h{k,x,t) = ±-[ x _f Q _ k a{\,t)e-^h{\x,t). (20) 

With the linear limit, obtained from the above solution by taking simply /i = 1, 
it can be verified that the evolution for q in (||) with 7 / has precisely the 
solution ([!(]). We remark at this stage that the time evolution of the nonlinear 
Fourier transform a(k, t) bears the linear character of the evolution. Hence 
the nonlinearity enters only in the output expressions (Q) and (^) through the 
nonlinear combinations of a with 7x,2- 

In order to realize how the mechanism of pump depletion and spike formation 
is allowed by the equation (|l6|), we deal with all quantities evaluated at k = 
(corresponding to a very sharp g(k)). Then we assume a Gaussian shaped laser 
pump input (from (||) the amplitude is normalized to 1) and a small proportion 
of Stokes seed with possibly a 7r-phase shift: 

h{t) = exp[-(i - hf/rf], I 2 {t) = tanh[(i - t ) / T ]h{t) /V . (21) 

With such input data, the pump repletion can be physically understood as 
resulting from a reversal of the Raman gain due to the change of sign of the 
Stokes input. This is precisely this behavior which is described by the formula 
( |l6| ) where, if h changes sign, then a starts to decrease and the pump repletion 
limit (|l^) is approaches. 

We have drawn the energy |ai| 2 of the pump output given by (0) on figs. 1 
and 2, for the above choices of h,i with the parameter values t\ = 50, t\ = 
45.5, T = 10, 3(0) = 100, 7 = 160 and no Stokes phase flip in fig 1 (that is 
to = 0) and a phase flip in to = 50 for the fig 2 (with ro = 1). These parameters 
correspond in the physical world to a pump pulse of maximum amplitude Aq of 
6.36 10 6 VmT 1 and a Stokes seed of 0.636 10 6 Vm^ 1 , for the values (taken from 
rcf. IH) for r in nanosec: N = 8.5 10 19 , \K LS \ = 5.8 lO" 17 and T 2 = 0.625 ns. 
Consequently, from 7 = L/icT^) = 160 the data would correspond to a beam 
path of 30 m. 
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PUMP DEPLETION 




40 60 
Time (ns) 



100 



Fig.l: Pump energy profile \a\\ 2 at x = +oo (input, dashed line) and at x = — oc 
(output, solid line) with unflipped phase of the Stokes seed. 



SPIKE OF PUMP RADIATION 
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Fig. 2: Pump energy profile \a\\ 2 at x = +oo (input, dashed line) and at x = — oc 
(output, solid line) when a ir-phase shift of the Stokes seed is introduced at 50ns. 

We can now easily understand the different behaviors of the pump output 
on figs. 1 and 2 just by inspection of (|l6|). Indeed, in the first zone (up to 30ns), 
a is small and we have no depletion. Then the growth of a (pump depletion) 
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results from the factor exp[0] in ( Jig ) because U given by (|ll]) is positive. Now 
if h in ( [TT] ) does not change sign (fig 1), a grows up to when the damping 
term dominates again and we observe the pump radiation again (right hand 
side hump in fig 1 and 2). If instead we introduce a phase flip in the Stokes 
seed, then the integral in (|ll]) makes a to decrease and possibly to vanish, 
which constitutes the mechanism for reversal of pump depletion. However this 
is allowed only if the growth of a with exp[</>] is not too fast, and here enters the 
role of the damping term exp[— "ft]. Hence the spike of pump radiation occurs 
as a result of Stokes phase flip via a balance between Stokes amplification and 
optical phonon damping. 

In conclusion, we have obtained the following set of results: 

1 - An exact and explicit solution to the SRS equations for short pulses (7 = 0). 

2 - An approximate solution for long pump pulses which describes perfectly in 
a unique formalism the pump depletion and the formation of a spike of pump 
radiation when the Stokes seed is given a phase flip. 

3 - The proof that the spike of pump radiation occurs as a balance between 
Raman gain and phonon damping, and that it can survive long propagation 
distances. 

4 - A new mathematical structure, the Raman spike, related to a zero of the 
reflection coefficient (our a(k,t)), and which, in the 1ST scheme, is part of the 
continuous spectrum. 

5 - The proof that the spike of pump radiation is not a soliton, namely that it 
is not related to a discrete part of the (nonlinear Fourier) spectrum. 

6 - An explicit formula for the description of transient SRS which can be used 
to study for instance the decay of the Raman spike, the generation of multi- 
spikes, the result of a stochastic phase in the generated Stokes wave, etc... These 
studies will be reported elsewhere (see [Q ) . 
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